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In this work the effect of the elbow-bend geometry and the effect of the tube arrangement on the 
performance of air-to-water heat exchanger is studied experimentally. In elbow-bend heat exchanger, 
the direction of the working fluid is bended at 90 degrees to its inlet direction. The heating or cooling 
fluid flows inside straight tubes while the working fluid flows past the tubes along an elbow pass. Three 
different types of the geometry of the elbow with three different tube bank arrangements were studied. 
The results were plotted and analyzed to clarify the effects of the elbow-bend geometry, the tube bank 
arrangements and the dead volume in the heat exchanger on the heat transfer and pressure drop. Two 
empirical correlations were deduced for each design, one to predict the relation between Nusselt and 
Reynolds numbers, while the other relation is between the friction factor and Reynolds number. This 
work was done to select the more suitable design to be used as a heater or a cooler in Stirling machines. 

Crown Copyright © 2010 Published by Elsevier Ltd. All rights reserved. 


1. Introduction 

Stirling engines can be operated using various heat sources. 
Referring to the literature of Stirling engines, different designs of 
heat exchangers are found. Some of them are explained below. A 
small gamma-type Stirling engine was tested by Iwamoto et al. [1 ]. 
In this engine the heat exchangers are located around the dis¬ 
placer cylinder. The heater and the cooler are of the fine-tube type. 
Hot water flows in the heater tubes. Coolant water flows in the 
cooler tubes. The working gas flows around these tubes. Hoshino 
[2] used twelve electrical heaters attached to a gamma-type Stir¬ 
ling engine heater head. The working gas is heated through 
annular flow passages of a groove-type heater; each of them has 
25 plate fins around a cylindrical heater. The cooler is of the shell- 
and-tube type, which has 322 tubes of 1 mm inner diameter. A 
heat input system which is composed of a heater and a burner as 
a newly devised heat exchanger in a Beta-type Stirling engine was 
designed and tested by Chang [3]. The heater likes a U-cup and it 
has slits which play the role of the heater tubes on the outer wall 
of the cup. A 200 W domestic free-piston Stirling electric power 
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system with multi-fuel capability was designed by Biao et al. [4]. 
The engine was fitted with a fine heater and an annular shaped 
clearance regenerator. The cylinder head was provided with 
vertical fins. At operation, the hot gases flow along the vertical 
fins, and the heat was supplied to the helium working gas through 
these fins. A thermodynamic analysis of a Stirling engine including 
dead volumes of hot space, cold space and regenerator was carried 
out by Kongtragool and Wongwises [5]. The analysis indicated that 
the engine net work is affected by only the dead volume while the 
efficiency is affected by both the regenerator effectiveness and 
dead volume. The heat transfer characteristics and pressure drop 
of block-type heat exchangers were experimentally investigated 
by El-Ehawany et al. [6]. The experimental data were compared to 
get the most suitable block-type heat exchanger to be used as 
a cooler and a heater in Stirling machines. Generalized correlations 
of heat transfer and friction factor for air-side plain fin-and-tube 
heat exchangers were reported by Wang et al. [7], where 31 
samples of fin-and-tube heat exchangers were used to develop the 
correlations. An experimental study to investigate the heat 
transfer and friction characteristics of fin-and-tube heat 
exchangers was carried out by Yan and Sheen. [8]. The results 
were presented as plots of friction factor and Colburn factor versus 
Reynolds number. Kim et al. [9] used a multiple regression tech¬ 
nique to correlate the data of 47 sets of plain fin-and-tube heat 
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Fig. 1 . a. Schematic diagram of the main parts of an alpha-type Stirling engine, b. Elbow-bend heat exchanger assembly, c. Elbow-bend heat exchanger, d. Isometric of the main parts 
of an alpha-type Stirling engine with elbow-bend heat exchangers. 


exchangers with staggered tube arrangements to develop the heat 
transfer and friction correlations. The geometry and the dimen¬ 
sions of the heat exchangers have a strong influence on the 
performance of Stirling engines. Therefore, the processes in the 
heat exchangers should be carefully considered in the modeling of 
the engine processes. Modern second-order thermodynamic 
model was used to simulate Stirling cycle to calculate the heat 
transfer and pressure drop during the engine cycle, and as 


a consequence, the obtained information is in turn used to rectify 
the prediction of the engine performance as provided by Mahka- 
mov and Ingham [10]. Moreover, recently, 2 and 3-dimensional 
CFD models have been developed for the modeling of the heat and 
mass transfer processes in heat exchangers of Stirling engines, 
which were considered as integral parts of the gas circuit as 
provided by Mahkamov [11,12]. The heat exchangers were also 
analyzed separately as given by Mahkamov and Ingham [13]. 
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1- Air inlet 

2- Air blower 

3- Air regulating valve 

4- Air duct 

5- Orifice plate 

6- Heating box 

7- Electric Heaters 

8- Mixing section 


9- Test section 

10- Test specimen 

11- Air outlet 

12- Water supply 

13- Water tank 

14- Water regulating valve 

15- Water inlet 

16- Water outlet 

Fig. 2. Experimental test rig. 


17- Electric power supply 

18- Voltage regulator 

19- Temperature measurement 
unit. 

20- Rotary selector switch 

21- Electric wires 

22- Stand 

23- To U-tube manometer 


The conventional designs of the alpha-type Stirling engine with 
twin parallel cylinders may have plain heater and cooler or tubular 
heater and cooler. For plain heat exchanger, insufficient heat transfer 
rates were observed and for tubular heat exchangers, the tubes are 
bent as a quadrant of a circle and the pressure drop is increased [14]. 
The suggested design aims at using straight tubes instead of bent 
ones. This restrains the flowing of the heat carrier fluid and the 
cooling agent through the tubes and the working fluid flows past the 
tubes as shown in Fig. 1. The major disadvantage of the present 
suggestion is that the dead volume of the engine could be increased 
compared with the traditional ones. However, the decrease in 
pressure drop is an important factor in the present aspect. 




Fig. 3. Isometric of the test specimen No. II. 


Fig. 4. Cross section of the test specimen inside the shell. 
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In the present work, an elbow-bend shell-and-tube heat 
exchanger was introduced as a new design for the heater and 
cooler of alpha-type Stirling engines. The elbow-bend heat 
exchanger consists of a shell in the form of a right-angle elbow 


and a tube bank fitted perpendicular to the flow in the elbow. 
Eight specimens of this heat exchanger were designed, manu¬ 
factured and tested experimentally. The experimental investiga¬ 
tion is revealed of the influence of the elbow-bend geometry, 
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Specimen No. (Ill), 136 tubes (circular Specimen No. (IV), 125 tubes (in-line 

arrangement), i// =45.2 % , S v = 345.44 m~' arrangement), i// =49.6 %, S v = 317.5 m~ l 

19 tubes in the outer row, S n = S p =7.5 mm 

S n =7.5 mm, S p =7.62 mm 


Fig. 6. Geometrical specifications of the test specimens No. I, II, III & IV. 
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Specimen No. (V), 142 tubes (in-line 
arrangement), y/ =47.9%, S v = 328.38 m~ ] 

Header: 2 x 20 x 100 + quadrant of 80 mm 
radius, S n =S p =7.5 mm 



Specimen No. (VII), 113 tubes (circular 
arrangement), y/ =54.4%, S v = 287.02 m 1 

18 tubes in the outer row, 

S n -8.18 mm, S p =8.51 mm 



Specimen No. (VI), 124 tubes (circular 
arrangement), y/ =50.0 % , S v - 314.96 m 

19 tubes in the outer row, 

S n —8.18 mm, S p =7.78 mm 


m 



Specimen No. (VIII), 111 tubes (circular 
arrangement), y/ =55.2 % , S v = 281.94m 1 

19 tubes in the outer row, 

S n =9.05 mm, S p =7.78 mm 


Fig. 7. Geometrical specifications of the test specimens No. V, VI, VII &VIII. 
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Fig. 8. Variation of the quantity (A p a x V d ) versus air mass flow rate. 


Fig. 9. Variation of the quantity (Qa/V d ) versus air mass flow rate. 
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tube bank arrangement, and dead volume (difference between 
the elbow volume and the volume occupied by the tubes) on the 
thermal performance and pressure drop of the elbow-bend heat 
exchanger. 

2. Experimental apparatus 

Fig. 2 shows a schematic diagram of the test rig used in this 
work. It consists of a wind tunnel system forcing hot air past the 
outside surface of the tube bank by a 3 HP centrifugal fan. The air 
duct was a circular duct of 52.3 mm or 104.6 mm inner diameter. 
The air flow rate was measured by using one of the two orifice 
plates (one for each duct) that were designed and manufactured 
according to BSI data catalog [15]. The air flow rate was controlled 
by a regulating gate valve. The air was heated by five heaters 
having a total power of 10.75 kW in a heating box. Four heaters 
were connected directly to the main supply, while the last one was 
connected through a voltage regulator. The air temperature was 
adjusted at 146 ± 1.5 °C at the inlet of the test section. The water 
was supplied through a tank with a float supported at 7 m above 
the test section to sustain a constant water flow rate during the 
experiments. The water flow rate was measured using a tank 
having a definite volume and a stop watch. The temperature was 
measured by thirty-one pre-calibrated T-type thermocouples. The 
inlet and exit temperatures of air were measured by two 


thermocouple meshes. Each mesh consists of 13 thermocouples 
connected in parallel to evaluate the average temperature [16]. 
The inlet and exit temperatures of the cooling water were 
measured by two thermocouples; one at inlet and the other at the 
water exit. The tube surface temperature was measured by using 
two thermocouples. The air temperature at the orifice plate was 
measured by one thermocouple. The thermocouples were con¬ 
nected to a digital temperature recorder model 922 (Manufac¬ 
turer: Maxtec Int. Co.). It has 0.1 °C resolution in a range up to 
200 °C. 

3. Test specimens 

Referring to Fig. lb, twin quadrants of a circular disc having 
5 mm thickness, were drilled using the plain radial drilling machine 
and the distribution of holes was accurately carried out by using the 
universal dividing head. A certain number of copper tubes having 
6.35 mm outer diameter, 0.4 mm wall thickness and 120 mm length 
were assembled to the twin quadrant by welding while keeping 
100 mm distance apart between the quadrants as shown in Fig. lc. 
Referring to Fig. 3, an isometric view of the specimen No. II is clear. 
Referring to Fig. 4, a square or a quadrant circle header was 
assembled to the twin discs. The tube arrangement allows the 
variation of both normal and parallel pitches s n ands p as shown in 
Fig. 5. In the vision of previous description, eight test specimens 
were designed and manufactured. The complete geometrical 
specifications of each one is given in Figs. 6 and 7. 
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Fig. 11. Variation of Nusselt number versus Reynolds number. 


Fig. 13. Variation of the quantity (Q^/V d ) versus air mass flow rate. 
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Fig. 14. Variation of the friction factor versus Reynolds number. 

4. Experimental data reduction 



Fig. 16. Comparison of present and previous studies. (NTU - s). 

NTU = U A/ C min (10) 


The analysis of the present work is concerned with the calcu¬ 
lation of non-dimensional numbers for each experiment and 
deducing generalized correlations to explain the case study. The 
following equations were used in the calculations: 
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The properties of the air and the water were found at the film 
temperature [17]. 
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5. Results and discussion 

In Stirling machines, increasing both or one of the dead volume 
and pressure drop of the heater or cooler causes deteriorations in 
the power and efficiency [7]. In the heat exchangers it is usually 
expected that the pressure drop increases as the dead volume 
decreases and visa versa. So, in the present work it is convenient to 
select the specimens that are having low (Ap a xl/ d ) (pressure drop 
times dead volume) and high (Qa/V d ) (heat transfer rate divided by 
dead volume). The quantities (Ap a x V d ) and (Qa/V d ) were 
compared for the specimens No. I, II, III, IV and V (the five speci¬ 
mens which have equal normal and parallel pitches) versus the air 
flow rate in Figs. 8 and 9. It is clear that the quantity (Ap a x V d ) 
tends to increase with the increase of the air mass flow rate for the 
five specimens. The specimen No. Ill exhibits the lowest value of 
(Ap a x V d ). The reduction in (Ap a x V d ) is due to its small dead 
volume V d , and it has a circular air pass about the tubes which 
minimizes eddies formation and as a consequence the pressure 
drop reduces. The specimen No. IV exhibits the highest value of 
(Qa/Vd). This is due to the streamlines of air not straight due to the 
circular pass of air through an in-line arrangement of tubes which 
provides different angles of attack between the air and the tubes. 
Referring to Fig. 10, the friction factor/is compared for the speci¬ 
mens No. I, II, III, IV and V versus Reynolds number. The specimen 
No. I exhibits the lowest friction factor/ due to minimum (Ap a ) for 
the same air flow rates. Fig. 11 implies that the specimen No. I gives 
higher Nusselt number compared with the other ones. 

From the pervious discussion, it is better to select the specimen 
No. Ill to be compared with the specimens No. VI, VII, and VIII. The 
quantities (Ap a x V d ) and (Qa/V d ) are compared for the specimens 
No. Ill, VI, VII, and VIII versus the air flow rate in Figs. 12 and 13. It is 
clear that the specimen No. VIII has the lowest (Ap a x V d ). Also, the 
specimens No. VI, VII, and VIII have approximately the same values 
of (Qa/Vd), that are slightly greater than that of specimen No. III. 
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Table 1 

Constants C\, C 2 , C 3> A and B. 


Specimen No. 

Range 

Ci 

c 2 

C 3 

A 

B 

I 

1100 < Re a < 2700 
2700 < Re a < 10500 
10500 < Re a < 10500 

0 

6.476 x 10 " 11 

0 

-3.174 x 10 ~ 5 
-1.646 x 10 ~ 6 
- 2.110 x 10~ 7 

0.1453 

0.06239 

0.05189 

9.327 x 10 " 3 

0.9389 

II 

1100 < Re a < 2700 
2700 < Re a < 8200 
8200 < Re a <11700 

0 

-1.480 x 10 10 

0 

-4.313 x 10 ~ 5 

9.567 x 10 ~ 8 

1.746 x 10 8 

0.1892 

0.07543 

0.06486 

3.691 x 10 " 3 

1.072 

III 

1090 < Re a < 2600 
2600 < Re a < 8900 
8900 < Re a < 12600 

0 

-1.379 x 10 10 

0 

-2.517 x 10 ~ 5 

1.297 x 10 ~ 7 
-3.680 x 10 -7 

0.1364 

0.06832 

0.05856 

3.598 x 10 " 2 

0.7995 

IV 

1100 < Re a < 2700 
2700 < Re a < 10600 
10600 < Re a < 12900 

0 

2.387 x lO " 10 

0 

-2.474 x 10 ~ 5 
-4.601 x 10 ~ 6 
-6.525 x 10 8 

0.1376 

0.07916 

0.05590 

2.203 x 10 " 2 

0.8520 

V 

1100 < Re a < 2600 
2600 < Re a < 9700 
9700 < Re a < 12600 

0 

1.881 x 10" 10 

0 

-2.525 x 10 ~ 5 
-4.238 x 10 6 
-3.207 x 10 " 8 

0.1377 

0.08029 

0.05596 

1.149 x 10 " 2 

0.9112 

VI 

800 < Re a < 2800 
2800 < Re a < 9700 
9700 < Re a <11900 

0 

-3.795 x 10 ' 10 

0 

-3.375 x 10 ~ 5 

2.546 x 10 " 6 
-2.348 x 10 ~ 7 

0.1580 

0.06765 

0.06107 

1.116 x 10" 2 

0.9478 

VII 

800 < Re a < 2800 
2800 < Re a < 8000 
8000 < Re a <11500 

0 

-1.018 x 10~ 9 

0 

-2.875 x 10 ~ 5 

7.084 x 10 " 6 
-3.628 x 10 ~ 7 

0.1552 

0.07237 

0.06937 

1.776 x 10 " 2 

0.9202 

VIII 

600 < Re a < 2200 
2200 < Re a < 7900 
8000 < Re a < 10,900 

0 

-6.458 x 10 10 

0 

-3.872 x 10 ~ 5 

5.183 x 10 ~ 6 
-5.862 x 10 ~ 7 

0.1289 

0.06403 

0.07171 

7.502 x 10 " 3 

1.069 

Referring to Fig. 14, the specimen No. Ill has the lowest friction 
factor /. The specimen No. VIII has the highest Nusselt number 

/ = C t x Re\ + C 2 

x Rf?a H - C 3 


(15) 


compared with the other specimens at the same Reynolds number 
and it shows adequate accordance with the diagram of the heat 


Nusselt Number can be correlated as: 


transfer as shown in Fig. 15. The effectiveness (e) was determined 
from equation (9) in which the heat transfer rate Qavg was deter¬ 
mined as the average of the air and water side values. Referring to 
Fig. 16, the effectiveness (e) was plotted versus the number of 
transfer units ( NTU ). This was done using the experimental results 
of all the heat exchangers which were tested in the present work. 


JVu = A x Ref (16) 

The constants Ci, C 2 , C 3 , A and B are shown in Table 1. 

6. Comparison of present and previous studies 


The relation was correlated as follows: 

£ = 0.6602 x NTU 0J527 (14) 

The friction factor can be correlated as: 


Fig. 16 shows a comparison between the values of the effec¬ 
tiveness (e) for all the heat exchangers in the present work and 
those calculated from the correlation e = 1 - e~ NTU in Ref. [17] 
assuming that C*~0. It can be seen that the experimental results 
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Fig. 17. Comparison of present and Jakob’s works. (Re - f). 
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Fig. 18. Flow pattern past an in-line-tube bank [18]. 



Fig. 19. a. Flow pattern past one tube [14]. b. Location of the next tube (white) in the present work. 


are very close to that in the previous work for NTU < 0.75. For 
higher values of NTU, the experimental results show higher values 
of effectiveness. 

Fig. 17 shows a comparison between the empirical correlation 
of the friction factor / which is given in Ref. [17] for in-line 
arrangements and the friction factor/of the specimens No. Ill, VI, 
VII, VIII. It is clear that the results of Ref. [17] have higher /than 



the present data. For an in-line tube bank, the location of the 
impact point depends on s p /d and Reynolds number. For s p /d = 6 
the point of attack and the point of impact coincide at cp = 0.0°, 
while for s p /d = 6, the point of impact is shifted to cp = 75°, as 
shown in Fig. 18 as provided by Kakac et al. [18]. 

The flow pattern past a tube is indicated by the image presented 
in Fig. 19a by Fraas and Ozisik [19]. Referring to this image, it is clear 
that, the eddy formation extended to a distance equals the tube 
diameter and hence the next tube should be located at a distance of 
S p < 1.6 d as shown in Fig. 19b. In Fig. 19b, the eddy formation will 
be dissipated and the point of impact is shifted to cp = 90° and as 
a consequence, there are no eddies on the rear surface of the first 
tube and no impact force on the front surface of the next one. The 
same can be applied on any two successive tubes. Then the inter¬ 
mediate surfaces of the tubes have no effect. The surfaces of the 
tubes in the longitudinal row are expected to be as a corrugated 
sheet as shown in Fig. 20. As a result, the pressure drop will be 
decreased. This is due to the values of s p /d are ranging from 1.20 to 
1.34 in the present work. This explains why the friction factor in the 
present work is less than that calculated from the correlation of 
Jakob, [17]. 

7. Conclusions 

The following conclusions were made from the experimental 
tests of the air-side elbow-bend geometry and water tube bank 
arrangements. 
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1. The circular arrangement of the tubes in elbow-bend heat 
exchangers is better than both in-line and staggered arrange¬ 
ments if both arrangements are having the same number of 
tubes. 

2. The number of tubes does not have significant effect on the 
pressure drop s p /d is less than 1.34. 

3. The tube bank arrangement through elbow-bend has a signifi¬ 
cant effect on the pressure drop. It increases if s n decreases and 
it also increases if s p increases. 

4. General correlations were deduced for the studied geometries 
of such heat exchangers that may be used in the design of 
alpha-type Stirling machines. 

Nomenclature 

A surface area, m 2 

B heat exchanger width, m 

C heat capacity rate, J/K 

C p specific heat at constant pressure, J/kg K 

d tube diameter, m 

f friction factor 

G mass flux of the air based on the minimum flow area, kg/s 

h heat transfer coefficient, W/m 2 K 

j Colburn factor 

k thermal conductivity, W/m I< 

L effective length of tube, m 

m mass flow rate, kg/s 

N number of tube 

Nu Nusselt number 

NTU number of transfer units 

p pressure, N/m 2 

Q. heat transfer rate, W 

Re Reynolds number 

S pitch, m 

S v surface to volume ratio, m _1 

T temperature, K 

U overall heat transfer coefficient, W/m 2 K 

V volume, m 3 

e heat exchanger effectiveness 

cp impact angle of stream lines, degree 

ji dynamic viscosity, kg m/s 

p density, kg/m 3 

\p porosity 

Subscripts 
a air 

avg average 

d dead 


i inlet conditions 

min minimum 

max maximum 

n normal 

o outlet conditions 

p parallel 

s surface 

t total 

w water 
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